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To describe the formation of clusters - in the general case a nonequilibrium process - 
in an interacting many-body system constitutes a new challenge for few-body methods. 
This may happen when the residual interaction between the quasi-particles leads to cor- 
relations. An example for such a system is nuclear matter. In the laboratory finite pieces 
of nuclear matter can be produced in heavy ion collisions. In astrophysics nuclear matter 
occurs e.g. during the supernova collapse and the formation of a neutron star. 

A microscopic approach to treat the formation of clusters uses a generalized quantum 
Boltzmann equation. This coupled equation has been numerically solved for nucleon f^, 
deuteron fy, triton f t , and helium-3 fh distributions utilizing the Boltzmann-Uehling- 
Uhlenbeck (BUU) approach PJ. The coupling between the different species is through the 
collision integrals /C[/jv, fd, ■ ■ ■)■ For the deuteron loss /Q ut (P, t), e.g., it is 

K. d ut (P,t) = Jd 3 kjd 3 k 1 d 3 k 2 d 3 h\(k 1 k 2 k 3 \U \kP)\ 2 dN ^ pnN 

xfN(h, t)f N (k 2 , t)f N (k 3 , t)f N (k, t) + ... (1) 

where Jn = (1 — /at)- The ellipsis denote further possible contributions, e.g. dd ^ tp, 
dd ^ hp or processes like yd # np, etc. The quantity U is the Nd — ► NNN break-up 
transition operator that in general depends on the medium. This dependence is neglected, 
if experimental cross sections are used to replace Uq in Eq. ([l]) — a standard technique 
and in many cases very successful. To calculate Uq including the self energy shift and the 
proper Pauli blocking and study the influence of the medium on different observablesa 
generalized Alt-Grassberger-Sandhas (AGS) equation has been derived earlier |3] [|. 
The effective few-body problem in matter arises in the Green function approach along 



with a cluster mean- field expansion ]T0| or the a Dyson equation approach [ |TT|j . 

Using the effective equations derived elsewhere f^-|8| we may study two important effects 
of the medium on the effective few-body-systems embedded in nuclear matter: (1) The 
change of binding energy, i.e. the self energy shift, (2) the changes in the reactions 
rates. Both effects are important and have consequences for the simulation of heavy ion 
collisions. 

The change of binding energy eventually leads to the Mott effect (where -Ebound — ► 0~). 
Not as dramatic as in Coulombic systems - here the Mott effect leads to the transition 
from isolating to conducting phase - it, however, influences the number of clusters and 
the energy spectrum produced in a heavy ion collision. The Mott density depends on 
the momentum of the cluster as for higher momenta blocking of the constituents of the 
cluster is less effective. The Mott momenta for deuteron and triton are shown in Fig. [l] for 
two different temperatures. Tritons are more stable than deuterons and at lower densities 
both clusters a more stable for higher temperatures. 
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Figure 1. Mott momenta for triton and 
deuteron. 
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Figure 2. Nd — > NNN break-up cross 
section for different densities. 



For a typical temperature of the heavy ion collision (in the final stage) we have cal- 
culated the in-medium cross section. The result is shown in Fig. |^. The threshold shift 
is because of the change of the deuteron's binding energy. A strong enhancement of the 
maximum of the cross section appears, however less change at higher energies. 

As a consequence the reaction time scales become much faster, when in-medium rates 
are used in the calculation instead of isolated ones. Fig. ||| shown the deuteron break-up 
time r bu (F, n N ) evaluated in linear response, where the life time of deuteron fluctuations 
5fd(P,t) depends on the deuteron momentum P and the nuclear density [|J. Similar, a 
chemical relaxation time T re \(nisf) can be defined, which results from a linearization of the 
respective rate equations [/j. The relaxation times are shown in Fig. [|. 

Finally considering a specific heavy ion collision it is possible to calculate the total 
number of deuterons coming out of a central collision of 129 Xe on 119 Sn at 50 MeV/ A 
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Figure 3. Deuteron break-up time as a 
function of the deuteron momentum. 



Figure 4. Chemical relaxation time as 
a function of nuclear density. 
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Figure 5. Total integrated number of 
deuterons function elapsed time. 
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Figure 6. Ratio of proton to deuteron 
numbers SIS db function of cm. energy. 



Fig. [5| shows the integrated number of deuterons taking into account gain and loss terms 
induced by the reactions in the medium. The net effect is a significant enhancement of 
the number of deuterons. Fig. || shows the influence of using in-medium rates on the 
spectrum of the proton to deuteron ratio along wuth experimental data. 

From the analysis is becomes clear that the medium modified elementary cross section 
and the proper self energy correction (binding energy shift, Mott effect) of the clusters 
should be included in the simulation of heavy ion collisions at that energies. 
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